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area-averaged mass-flow coefficient on porous model surface 
cross section area of porosity probe tube 

area-averaged function of local mass-flow coefficient and surface 
pressure on porous model surface 


stagnation point mass-flow heating rate parameter, — 


mAH 


qo 


mass-flow heating rate parameter, reference 9 


mAH 


area-averaged mass-flow heating rate parameter, — 

^os 

local mass-flow coefficient on porous model surface 

skin friction coefficient 

specific heat of air or calorimeter slug 

roughness height 

local mass-flow addition rate 

free-stream mass-flow rate 

area-averaged mass-flow addition rate to a given location, s 

• 

m 

mass-flow rate ratio, -s — 

ni_ 


pressure 

reference heat transfer coefficient (stagnation point value on 
hemisphere model at m = 0) 

local heat transfer coefficient at s with m = 0 


area-averaged heat transfer coefficient up to location s on model 
surface with m = 0 

local heat transfer coefficient at s with m ^ 0 


Re Reynolds number 

r normal distance to point on model surface from axis of symmetry 
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s radial distance from model apex (stagnation point) measured along 

model surface 

t time, sec 

T temperature, K 

T’ turbulence intensity, reference 8 

AH heat-transfer driving potential, - T^) 

0 momentum thickness 

y viscosity of air 

p air density or calorimeter density 

T calorimeter slug length 

^ azimuthal angle on model surface 

^ stagnation point heating rate ratio, 


function of B’ , reference 9 

‘Is 

heating rate ratio, — 


Subscripts : 

D model diameter 

ex external surface of model 

in internal surface of model 

lam laminar 

St stagnation point 

trns transition point on model surface 

tt tunnel total condition 

turb turbulent 

w model wall or surface 

oo free-stream condition 
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EFFECTS OF MASS ADDITION ON BLUNT-BODY BOUNDARY-LAYER 


TRANSITION AND HEAT TRANSFER 
George E, Kaattari 
Ames Research Center 


SUMMARY 


Heat-transfer data were obtained on blunt models at Mach number 7.32 and 
free-stream Reynolds numbers in the range of 0.6x10^ to 5.2x10^ to investi- 
gate the effect of ablation on boundary-layer transition. Ablation was simu- 
lated by mass addition of air through the porous surface of the models. The 
ratio of the mass addition rate to the free-stream mass flow was varied from 
0 to 0.5. The models were 17.8 cm in diameter and consisted of hemispheres, 
spherical segments, and blunted 60° cones. 

The data were compared with various applicable boundary-layer codes in 
the laminar and transitional flow regimes. Empirical heating rate data cor- 
relations were developed for the laminar and turbulent flow regimes. 


INTRODUCTION 


Vehicles entering planetary atmospheres at high speeds are enveloped by 
high-temperature shock layers. At sufficiently high speeds, certain surface 
materials may ablate from the vehicle in gaseous form and protect the vehicle 
at the expense of surface material mass loss. The transferred heat from the 
high-temperature shock envelope is expended by converting a thin layer of the 
vehicle surface into hot ablative products that are subsequently carried away 
by the external flow field. Problems then arise in predicting the effect of 
ablation of heat transfer (particularly in the turbulent flow regime) and in 
predicting the effect of ablation on the transition from laminar to turbulent 
flow. 


Much experimental and theoretical study has been devoted to predicting 
laminar heat transfer in the presence of ablation for both blunt and slender 
configurations. These studies (refs. 1-5) indicate a reduction in heat trans- 
fer and have produced correlations and theory relating the heat transfer 
decrease as a function of the ablation rate. Along with ablation, the effects 
of surface roughness and free-stream turbulence on transitional and turbulent 
flow have been studied (refs. 6-9). Some results of these investigations 
are compared with the results of the present investigation. 

Studies of heat transfer in the literature are largely restricted to low 
ablation rates. The present investigation involves massive rates, approaching 
those anticipated for Jovian atmosphere entry. The primary purpose, along 



with providing boundary-layer transition data, is to determine to what extent 
massive ablation affects turbulent heat transfer rates. This aspect of abla- 
tion is important both in planetary exploration and in military application 
since, if the heating rate is significantly increased, nonuniform heat shield 
erosion and attendant thermal stress and aerodynamic stability problems may 
arise. 


EXPERIMENTAL CONSIDERATIONS 


Models 

The models used consisted of three 6. 35-mm-thick porous, stainless steel 
headers in the forms of a hemisphere, a 60° blunt cone, and a spherical seg- 
ment, as shown in figure 1. The headers were made interchangeable to a steel 
pressure chamber to which air was led through piping. The assembled header, 
chamber, and support system are shown in figure 2. Each header was provided 
with 15 to 16 calorimeters. A calorimeter (fig- 1) consisted of a 1.6-mm- 
diameter, 1 . 27-Tnm-long copper slug tightly fitted into a cavity at the end of 
a 2.4-iran-diameter, 13.7-mm-long insulating dowel. A chromel-constantan junc- 
tion was peened to the base of the slug and the wire leads passed through the 
dowel to a male plug. Each calorimeter was press-fitted to the header hole so 
that the dowel and slug surface were flush to the header external surface. 

Each header was also provided with five 1.6-mm stainless steel tube pressure 
taps. Figure 1 gives the location of the calorimeters and pressure taps for 
each header. Although the figure indicates spacing of the calorimeters in a 
single radial plane, they were actually arranged in a spiral 360° pattern to 
reduce "downwind” interference effects. Not shown in figure 1 are thermo- 
couples that were interdigitated with the calorimeter locations and spotted to 
the internal surface of the header. These thermocouples monitored the enter- 
ing air temperature during the tests. The pressure taps were located in a 
single radial plane as indicated. 

The models were intended to have uniform porosity. No attempt was made 
to "tailor" header porosity to simulate ablated mass addition rate distribu- 
tion of actual atmospheric entry. Pretest porosity surveys of the headers 
indicated the presence of manufacturing imperfections resulting in deviation 
from uniform porosity. Figure 3 shows a plot of representative mass-flow rate 
variations with radial distance s determined from the porosity surveys. 

Except for deficiencies near the stagnation region existing for the hemispheri- 
cal models and the conical models (not shown), the porosity deviation was 
generally within ±10% of the mean value. 

Initial tests were made to check out the systems and to determine the 
feasibility of the thermocouple installation. Some model degradation occurred 
during these initial tests. The hemispherical and spherical segment models 
were refurbished, and the conical model was replaced since the original model 
proved to have a nearly impervious nose section. The refurbishing and 
replacement resulted in models with somewhat different porosity distributions. 
The altered models were tested to a higher mass^flow rate range in the second 
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phase of the investigation. To distinguish the results presented herein, the 
models are referred to by phase number, that is, by model 1 or model 2. 


Test Conditions and Procedures 

The tests were performed in the Ames 3.5-Foot Hypersonic Wind Tunnel at 
Mach number 7.32. A description of the facility and its operations appears 
in reference 11. The Reynolds number, based on model diameter, ranged from 
0.6x10^ to 5.2x10^ and the temperature from 660 to 880 K. 

A test run established a predetermined flow rate of air through the 
model at about 27° C, the flow being gaged with a swirlmeter. The model was 
then inserted into the tunnel winds tream for about 1-1/2 sec, during which 
interval a shadowgraph picture of the model was taken and the time-temperature 
history of the model calorimeters was recorded. The pressure transducers 
reached equilibrium well within the insertion time. The mass-flow rate ratios 
with respect to free-s tream flow varied from 0.0 to a maximum value of about 
0.5 for cone model 2 at the lower P^eynolds number range. 


Data Reduction 

Mass-addition rate distribution— Pretest porosity surveys of the model 
headers were required to determine mass-addition flow rates as a function of 
radial location, s. This was accomplished by pressurizing the chamber and 
holding to the external header surface a 6.4-inm-diameter ^sniffer” tube 
attached to a low-capacity floatmeter. Floatmeter readings were taken at 
6.4-mm intervals in the radial direction and at 5° intervals in the azimuthal 
direction. The local conductance Cg at each location on the header surface 
was then evaluated with the equation 

m^ - CgA(P^^ - P|x^ 

where 

m = flow rate measured by the floatmeter 
A = area of the ^sniffer” tube 
P^^ = chamber pressure 

P^^ = external (atmosphere) pressure 


The local conductances were area-integrated to determine the mean conduc- 
tance. This conductance coefficient was then checked with a high-capacity 
venturi meter by passing air through the header. The resulting flow rate was 
compared with that calculated by equation (1) using the mean conductance and 
total header area. The comparison gave satisfactory agreement for a wide 
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range of chamber pressures. Proportional adjustments to the local values of 
Cg were made to bring the calculated flow rate into exact agreement with the 
venturi-meter measurements. 

The mass-addition flow rate distribution used in this investigation was 
defined as the area-integrated unit mass flow to the radial position s of 
interest. The mass flow rate determinations for the test runs involved a 
variable external pressure as well as the variable conductance coeffi- 

cient c . These were taken into account by equation (2). 
s 


m = a P? 
— s in 


b P^ 
s tt 


( 2 ) 


where P^^ is the tunnel total pressure. The coefficients a^ and b^ 
were tabulated as functions of radial location s, dependent on the local 
conductance of each model and the external pressure suitably related to the 
tunnel total pressure. Detailed derivation of equation (2) is given in the 
appendix. 

The mass-flow rate distributions for all test runs of the investigation, 
calculated with equation (2), are presented in tables 1 to 4, 

Heating rate distrihutiovr- The heat transfer rates were calculated from 
calorimeter temperature versus time transients recorded while the model was in 
the tunnel stream. The usual formula was applied, namely 


= PV (d^) 

Because of slight misalignments in positioning the calorimeter surfaces 
flush to the external header surface, repeatable scatter in the heat trans- 
fer rates was noted. The magnitude of the scatter was determined by comparing 
the no-flow results of the present tests with unpublished test results of 
identically proportioned, solid thin-skin models tested in the same tunnel and 
flow conditions. Correction factors were determined for each calorimeter such 
that the ablation model results at the condition of no-flow were made to agree 
with those of the solid thin-skin models. 

Transition location— The tested models had a decreasing heating rate dis- 
tribution with distance s from the stagnation point at the lower Reynolds 
numbers and at near-zero mass addition rates. This distribution is character- 
istic of laminar heating for the model shapes considered. The transition from 
laminar to turbulent flow at certain increased mass addition rates was made 
evident by sudden increases in heating rates occurring beyond a definite loca- 
tion s. This location or transition point in the present investiga- 

tion was defined as the last location beyond which point a reversal to an 
increased heating rate occurred. The transition point is evident in fig- 
ures 4(a)- (e). The transition point and other significant features of the 
heating rate distributions are labeled in figure 4(c) where they are 
particularly well defined. 
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RESULTS AND DISCUSSION 


Experiment 

The heat-transfer distribution results for the range of variables test 
are presented in figures 4 to 9 in the form of normalized heat-transfer 
coefficients plotted as a function of distance s from the stagnation point 
location of the models. A nominal mass addition rate ratio parameter m* 
is associated with each curve of the figures. The test numbers and pertinent 
flow variables are also indicated. The reference heat-transfer coefficient 
is the theoretical stagnation point value for a 17 .8-cm-diameter hemi- 
sphere with no mass addition at the test conditions of each run. 

Table 1 contains tabulations of the mass addition flow rate distributions 
m as a function of s for each run presented in figures 4 to 9. The value 
of the free-stream mass-flow rate ra^o of each run is also listed for con- 
venience in determining the local mass addition flow rate ratio (m* = m/m^) 
if desired. 

Hemispherical models— The results for hemispherical model 1 are presented 
in figure 4. The mass addition initially decreased the level of heating over 
the entire model surface with increasing m and indicated the presence of a 
laminar boundary layer. Then, with further increases in a transition to 

a turbulent boundary layer occurred, as indicated by a rapid increase in heat- 
i ing rates beginning at definite locations s as previously described 

(at s = 2.5 cm in fig. 4(a) and at s = 5.8 cm in fig. 4(c)). The transi- 
tion location shifted progressively to lower values of s with increasing 
values of m*, as is particularly evident from figures 4(c) and 4(e). Con- 
tinued increases in m* ultimately decreased the heating rates, as is shown 
in figure 4(d). Note that the heating rates at, and near, the stagnation 
point s = 0 are relatively insensitive to m*. This phenomenon was also 
noted in a similar investigation of a hemisphere reported in reference 1. 

The heating rate distributions for hemispherical model 2 are presented 
in figure 5. Most of the data were taken at higher mass addition flow rate 
ratios m* than with model 1. Tests with model 2 were also run at a lower 
Reynolds number (0.5x10^). The mass addition rate distribution of the two 
models (tables 1 and 2) differs somewhat, model 2 being less permeable to 
flow near the stagnation point. Model 2 also had a local flow rate irregular- 
ity at s ^ 1.3 cm, which caused a conspicuous perturbation in the heating 
rate in that vicinity. The transition point movement with increasing m^^^ is 
apparent in figure 5(e). 

Spherical segment models- The results for the spherical segment models 
are presented in figure 6. The effect of small negative mass addition rates, 
shown in the data of figures 6(a) and 6(c), causes large increases in heating 
rates, undoubtedly due to partial inspiration of the oncoming hot airstream. 
Positive values of mass addition affected the heat transfer in a similar 
manner as that noted for the hemispherical models. In the case of the spheri- 
cal model, however, no progressive movement toward the stagnation region 
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occurred at the lowest Reynolds numbers (fig. 6(a)) and the transition loca- 
tion became established at s = 3.3 cm and remained fixed for values of m* 
from 0.018 to 0.079. A perceptible shift of the transition point from 
s = 3.6 cm to s = 2.5 cm occurs in the range of m* presented in the 
intermediate Reynolds number range (fig. 6(b)). The largest Reynolds number 
data (fig. 6(c)) exhibit the same behavior as those of figure 6b; however, 
calorimeter malfunction precluded definite fixes for the transition location. 
The large increase in heating rate in the stagnation regions at Reynolds 
number 4.93x10^ with m* = 0.004 (round symbols, fig. 6(c)) was probably 
due to combined effects of mass addition and tunnel free-stream turbulence 
since, at this Reynolds number, enhanced heating rates were also found at 
m* = 0 in preliminary tests of a thin-skin nonporous version of the spherical 
segment model (dashed-line curve) . The spherical segment model did not 
exhibit the stagnation point heating rate insensitivity to mass addition to 
the degree noted for the hemispherical models. 

The heating rate distributions for the spherical segment model 2 are pre- 
sented in figure 7. Model 2 generally exhibited higher heating rates near 
the stagnation point than did model 1 under similar flow conditions. The 
difference in heating rates was probably due to differences in the mass addi- 
tion rate distributions between the models (tables 2 and 4) . Model 2 had 
a somewhat irregular mass flow distribution, with minimum values near s = 0. 
Model 1 had a nearly constant mass flow distribution over the entire model 
surface. The nonuniform mass addition rate distribution of model 2 may have 
caused its somewhat earlier transition to turbulent flow. 

Coniodl models— The test results for conical model 1 are presented in 
figure 8. The effects of mass addition on the heating rate distribution fol- 
lowed the pattern of the hemispherical models. The transition point moved 
forward (towards s = 0) with increasing values of m’*^ (figs. 8(a), 8(c), 
and 8(e)). Large heating rate increases resulted when transition-to-turbulent 
flow was attained at the larger Reynolds numbers (fig. 8(c), run 63C1 and 
fig. 8(e), run 70C1). The range of s which was relatively impervious to 
flow is noted (0 < s < 3.4) in figure 8. However, there is a small effect of 
mass "leakage” in this region which, at the highest Reynolds number (fig 8(e)), 
was sufficient to cause early transition to turbulent flow (runs 71C1 and 
73C1). 

The heating rate distributions for cone model 2 are presented in 
figure 9. Comparison of the data between models 1 and 2 at approximately the 
same flow conditions (run 45C1 and 42C2) indicates turbulent flow for model 2 
and transitional flow for model 1. The differences in heat transfer rates 
between the models must then be due to gross differences in mass addition rate 
distributions (tables 5 and 6). Model 1 had an nearly impervious nose sec- 
tion; model 2 was porous over the entire surface. Table I indicates, however, 
that the mass flow rate distribution for model 2 varied in an irregular manner. 
A particularly "high spot" occurred at s = 1.27 cm. This is reflected in the 
heating rate data of figure 9 by a substantially reduced heating rate at this 
location. It should be clarified that the irregular mass flow distribution 
of model 2 was in respect to azimuthal location. The radial or s-wise distri- 
bution was fairly uniform upstream of each calorimeter location; thus, "high 
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spots” of mass addition were only apparent locally and were not instrumental 
in promoting early flow transition. 


Correlations 


The results of the investigation were examined and attempts were made to 
correlate the results. Limited comparisons were made with the results of 
other investigations. 

Lconinar heating rates-- A semi empirical method for correlating the effect 
of mass addition on laminar heat transfer distribution is described and com- 
pared with a more sophisticated method and with experimental results of the 
present investigation . 


The effect of mass addition on stagnation point heating (ref. 2) is 
generally well predicted for laminar flow by 



= 1 - 0.72B + 0.13b2 


(4) 


The term ip is the ratio of the heat transfer coefficient qg at the 

stagnation point in the presence of mass addition with respect to the heat 

transfer coefficient with no mass addition. The independent variable 

B = mAH/q^ is the mass addition parameter. Equation (4) was found to apply 

to locations other than the stagnation point when the term ip was redefined 

as ^ = q /q , or the ratio of the heat transfer coefficient at location s 
s os 

in the presence of mass addition with respect to the coefficient at the same 
location with no mass addition. The mass addition parameter B is modified 
to ^ = MH/qos* The term m is the area-averaged mass-addition rate to the 
location s and 1/qos Ts the area-averaged, inverse heat transfer coeffi- 
cient (no mass addition) to the same location 


q 


-1 

os 



(5) 


Equation (4) in the modified form was tested for the hemispherical model 
by correlating surface heating rates generated with the boundary-layer code 
reported in reference 5. The results are shown in figure 10(a) for various 
mass addition rate distributions. The agreement is considered good. Some 
small second-order deviations from perfect agreement are apparent. 


Experimentally determined surface heating data for the hemispherical 
models were similarly correlated with equation (4). The results are shown 
in figure 10(b). Note that the conspicuous perturbations in the heat trans- 
fer coefficients of model 2 at s ^ 1.7 cm (shown previously in fig. 5) do 
not impair the correlation of figure 10(b) (solid symbols). 
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The correlation of the heating data over the spherical segment model 
surfaces with modified equation (4) is shown in figure 10(c). Generally fair 
agreement is indicated except for the data of run 97SS2. 

The data of conical model 1 (fig. 10(d)) correlates well, with the 
exception of some data points of run 45C1 and an isolated data point of run 
63C1. Only the data of model 1 were considered since most of the model 2 data 
were taken at high mass-addition rates where transition occurred close to the 
stagnation point. 

No effect of Reynolds number was apparent in the correlations of the 
data of figures 10. 

Anomalous stagnation point heating— In the experimental results, the 
heating rates at the stagnation point were relatively insensitive to mass 
addition rates. This fact was most evident for the hemispherical model data. 
In reference I, similar results were noted for a hemispherical model tested 
at Mach number 5 with nitrogen gas for mass addition. The results of refer- 
ence 1 and the present data are compared in figure 11(a). The free-stream 
Reynolds number based on model diameter is closely matched. Statistically, 
the present results indicate less effect of mass addition than those of refer- 
ence 1, although the results overlap at Reynolds number 5.2x10^. No con- 
sistent effect of free-stream Reynolds number is evident in figure 11(a) in 
the data of reference 1 or of the present investigation. 

The effect of mass addition on stagnation point heating for the spherical 
segment model is shown in figure 11(b). The data indicate that the heat 
transfer rates are reduced significantly compared to those for the hemispheri- 
cal model. The data tend to approach the laminar distribution, but neverthe- 
less are still far above theoretical expectations at large injection rates. 

The effect of mass addition on the stagnation point heating of the coni- 
cal models is not well documented. Model 1 was impervious to mass addition 

at the stagnation point and only three test runs performed in the mass addi- 

tion range, 0 < ]B < 3, with conical model 2. These runs gave values of 
q/q^ ranging from 1.01 to 0.95, indicating the least effect of mass addition 
on the stagnation point heating among the three configurations tested. 

The apparent configuration effect on the stagnation point heating indi- 

cated by the present data is probably a scale effect due to nose radius, 
which primarily affects the stagnation point velocity gradient. 

High heating rates (above laminar) at the stagnation point can be caused 
by free-stream turbulence (ref. 6). The difference in the two sets of data 
in figure 11(a) may then be ascribed to the turbulence levels in the facili- 
ties used. A peculiarity of the heating rate distribution immediately down- 
stream of the stagnation point occurs in the present data and is particularly 
evident for the hemispherical models in the ^ coordinate plots of 

figure 12. The data of the indicated runs do not correlate with the pre- 
dicted laminar (dash line) curve, but plot on a series of parallel curves. 
Apparently, the enhancement due to turbulence at the stagnation point does 
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not persist downstream, and the flow tends to be laminar-like up to the point 
where transition ultimately occurs. 

Boundary- layer transition— Gaseous mass addition from a body surface 
forms a layer that interfaces with the external flow. At small mass addition 
rates, the layer is thin and the interface is stable. When mass addition is 
increased, the layer thickens and, ultimately, the interface becomes unstable. 
Turbulent mixing with the external flow ensues with concomitant increase in 
heat transfer to the body surface. The beginning of instability, or transi- 
tion, was evident in the data previously described where a sudden increase in 
heating rates occurred with increased mass addition. The problem of predict- 
ing the occurrence of transition as affected by mass addition is complicated 
by other contributing factors. Two significant parameters have received 
experimental and theoretical study (refs. 6-9) and are discussed in the light 
of the present investigation. 

Effect of free-stream turbulence: While free-stream turbulence was not 

an experimental variable in this investigation, it is an important factor in 
comparing data from different facilities. A correlation between computed 
transition Reynolds numbers utilizing the test conditions and tunnel turbu- 
lence level T of the present investigation was found by Mr. Wilcox of DCW 
Industries and is presented in reference 8. This theoretically based corre- 
lation relates the boundary-layer edge Reynolds number to the local flow rate 
momentum thickness Reynolds number at the transition location of the hemi- 
spherical models by 


trns 

43000 


1+41 exp 


/_ 13 ^^trns\ 
\ 12 yw y 


( 6 ) 


Experimental transition point Reynolds numbers for the hemispherical 
models are shown in figure 13. The data are represented by sets of two points. 
The lower point of each set indicates incipient transition. The upper point 
(with arrowhead) represents completed transition. The data sets bracket the 
theoretical correlation curve at all Reynolds numbers. However, in general, 
the experimental values indicate earlier transition than is predicted, 
particularly at the lower Reynolds numbers. The relatively poor correlation 
at the lower Reynolds numbers (lower surface pressures) may possibly be 
attributed to ”mass-flow roughness” or minute jetting action, which tends to 
be suppressed at the higher Reynolds numbers. 

Effect of surface roughness: A degree of surface roughness existed in 

the models of the present investigation because of the permeable, porous 
structure necessary to permit mass flow through the surface. A photomicro- 
graph of the hemispherical model surface section (fig. 14) shows a fairly 
consistent surface texture in terms of ”valley-to-peak” height. This dimen- 
sion was taken as the roughness parameter k of the model. The experimental 
transition Reynolds number based on momentum thickness was compared with those 
given by the correlation relationship of reference 9, where combined roughness 
and mass addition rates are taken into account 
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where 


Qualifications for application of the above correlation relationship are 
stated in reference 9: ”A computed value of 255 must be reached or exceeded 

at the sonic point location; if this condition is satisfied, the transition 
zone is predicted to physically begin at the surface point where this param- 
eter attains a value of 215. If predicted to occur, the transition zone will 
always be located in the subsonic flow region." The correlation relationship 
(eq. (7)) was tested with the hemispherical model data. The results are 
shown in figure 15. The data group labeled Re 0 < 255 and two data points 
with transition occurring beyond the sonic point do not meet the qualifica- 
tions quoted above and also do not correlate. The remaining data meet the 
qualifications and are in fair accord with the indicated PANT correlation 
curve . 

Turbulent heating rates— The ^ curves of figures 10(a) through 

10 (d), discussed previously, relate the heating rate in the presence of mass 
addition with reference to a no-mass-addition laminar heating rate as a 
function of the mass flow rate parameter. Satisfactory correlation of the 
experimental results over the model surfaces was obtained at low mass-addition 
rates. Logically, the effect of mass-addition in the cases where turbulent 
heating is predominant should be measured with respect to a no-mass-addition 
turbulent heating rate reference. A correlation of the data on this premise 
was made and is described in the following paragraph. 

Turbulent flow factor: The procedure for correlating the turbulent 

heating rate data of this investigation was simply to multiply the laminar 
based coordinates and B by the factor qo/^oturb’ thus shifting them to 
a turbulent heating rate reference. Numerical values for the factor were 
obtained by applying Reynolds analogy (proportionality) between heat and skin 
friction coefficients 



^lam ^ 1.328xRe ^ 

0.072xRe“°‘^ 

turb 

The numerical values for the skin friction 
reference 10 . 


18.44 


Re 


0 . 3 


""turb 

coefficients Cf 


( 8 ) 


were taken from 


Turbulent heat transfer correlations: The turbulent flow factor defined 

in equation ( 8 ) was applied to the laminar parameters ij; and ^ to convert 
them to turbulent flow parameters 
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(i^), 

<x> 

D 

and 



D 


Turbulent data were taken as the values downstream of the peak heating rate 
after transition (see fig. 4(c)). The turbulent data of the hemispherical 
models were plotted in the above defined J^turb " — turb coordinates of 
figure 16(a) as case I. All data of both models and Reynolds numbers corre- 
late except for the data of model 1 at the highest Reynolds number (solid 
symbols) . The data from the hemispherical model at this Reynolds number 
(fig. 4(e)) show a more pronounced shift of the transition point with mass 
addition than is generally the case with the data at lower Reynolds numbers. 

An attempt to account for the fact that the flow was not fully turbulent 
over the entire model surface in all cases was made in the following approxi- 
mate manner: The local reference q^g (eq. (5)) was assumed to retain 

laminar values up to the midpoint location between the transition point and 
peak turbulent heating point and then to assume the fully turbulent value 
to the point in question. The ’’weighted" reference ^Qg was used to modify 
the values of ^turb ^turb* data, plotted in the modified coordi- 

nates, are shown as case II in figure 16(a). The high Reynolds number data 
in this case correlates with the data at the lower values. The scatter of the 
data is random and the correlation is considered satisfactory. 

The data of conical model 1 was plotted only in the fully turbulent 
coordinates defined by equation (9). The results are sho\^^l in figure 16(b). 

An excellent correlation of the data in the Reynolds number range tested is 
evident. The high degree of correlation is due, in part, to the fact that 
the plotted data were confined to runs having approximately the same transi- 
tion point location. In addition, the reference coefficients qQg and qQg 
are relatively independent of s in the transition region; thus, ^turb 
^turb less sensitive to transition point shifts than in the c"ase of the 

hemispherical model. In this regard, preliminary plots of the data in 
"transition modified" coordinates indicated no significant changes to the 
degree of correlation shown by the fully turbulent correlation shown in 
figure 16(b). 

Correlation of conical model 2 and the spherical model data was not 
attempted. Thermocouple failures in the region of anticipated turbulent flow 
for conical model 2 precluded well-defined indications of turbulent heating. 
Data for the spherical models did not indicate the development of fully 
turbulent flow as was well-noted by the "peak heating" locations evident 
for the hemispherical models and conical model 1. 


lb 

—turb 


( 9 ) 
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CONCLUSIONS 


An experimental investigation was made of the effects of mass addition on 
the heat transfer to hemispheres, spherical segments, and blunted 60° cone 
models at Mach number of 7.4 in free-stream Reynolds number range 0.6 to 
5.2x10^ based on model diameter. The ratio of mass-addition rate to free- 
stream mass flow was varied from 0 to about 0.5. The results of this investi- 
gation were compared with empirically derived correlations and with the appli- 
cable results of various computer codes. Significant features of the data and 
conclusions drawn from comparisons of the data with applicable boundary-layer 
theories follow. 

The theoretically predicted effect of mass addition on stagnation point 
heating in laminar flow (ref. 2) was extended by an empirical, area-averaged 
heating reference technique to correlate laminar heat transfer distributions 
over the models* surfaces. The method was in good agreement with the results 
of both experiment and the laminar boundary-layer code of reference 5. 

Anomalously high stagnation point heating based on laminar theory was 
noted in the data of the present investigation. The same results were found 
in a similar investigation reported in reference 1. Enhanced stagnation point 
heating is generally believed to be due to free-stream turbulence. The data 
of the present investigation at low mass addition rates seem to indicate that 
the anomalous stagnation point heating is due to localized turbulence with 
subsequent downstream flow tending to become laminar. 

The effect of free-stream turbulence on boundary-layer transition has 
been investigated theoretically (ref. 7). Using this work as a basis, and 
utilizing the tunnel test conditions and turbulence level of the present 
investigation, reference 8 presents a formula (eq. (6)) relating the boundary- 
layer edge transition Reynolds number to the local momentum thickness and 
ablation rate for the hemispherical model. This relationship correlated the 
experimental data satisfactorily, particularly at the higher Reynolds number 
where transition was markedly abrupt. 

The effect of surface roughness on transition has been investigated 
theoretically (ref. 9). This work presents a correlation of momentum thick- 
ness Reynolds number at transition as a function of local mass flow, momentum 
thickness, and roughness parameter (eq. (7)). Although not a considered 
variable in the present investigation, roughness measurements were available 
for the hemispherical model and some experimental confirmation of the corre- 
lation was found. The correlation was limited in extent due to qualifying 
restrictions . 
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The empirical correlation parameters ^ and B for laminar flow heat- 
transfer distributions were modified by application of Reynolds analogy, and 
correlation of heat-transfer distribution data was extended to areas of 
turbulent flow on the conical and hemispherical models. 

Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, California 94035, Sept. 8, 1977 



APPENDIX 


DERIVATION OF MASS-FLOW RATES 


The computation required for the area-averaged mass-flow rate distribu- 
tion for each test run was cast into a convenient form requiring only the 
measured model internal pressure, the tunnel total pressure, and tabulated 
factors predetermined for each model. The derivation of the form used in the 
computation follows. 


The local mass-flow rate associated with the differential area element, 
rd(})ds, at the radial distance, r, from the stagnation point is 


m = 


c (Pf 
s xn 


P^ )r d(j) ds 
ex 


r d(|) ds 


(Al) 


Pj_j^ and P^^ are the local internal and external pressures on the porous 
header, and Cg is the flow coefficient at location r. 

Equation (Al) is integrated in the azimuthal direction . 


d(() 


m = 


/: 




d(|) 


/; 


(A2) 


r ds 


The value of P^^^ is constant for a given run. The variable P^^ 
a function of s and of the tunnel total pressure 


IS 


P = f(s,M )P^^ 
ex ’ oo tt 


from 


since 


and 


P ^st 


ex 


St tt 


p 

P,-. tt 


= f(s) , theory or experiment 


St 


St 


tt 


= f(M ) , a constant 

^ cry ^ ^ 


Substituting the above into equation (A2) 

's 


m = 


/ c [P? - f^(s,M )P^. ]r ds 

lo 3 p2 _ b p2 

s , s in s tt 


/: 


(A3) 


r ds 
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TABLE 2.- MASS ADDITION RATE DISTRIBUTION, HEMISPHERICAL MODEL 
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TABLE 4.- MASS ADDITION RATE DISTRIBUTION, SPHERICAL SEGMENT MODEL 2 
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Figure 3.- Typical variation of mass flow rate with radial distance s with 

constant pressure drop across headers. 
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Figure 4.- Concluded. 
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Figure 5.- Continued. 
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Figure 5.- Continued. 
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Figure 6.- Continued. 
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s, cm 


(a) = 138 N/cm2 (200 psi) , -0.001 < m* < 0.487. 

Figure 7.- Heat transfer distribution over spherical segment 
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Figure 7.- Continued. 
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(d) Pj.j. = 552 N/cm2 (800 psi) , -0.003 < m* < 0.020. 
Figure 7.- Continued. 
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(a) = 276 N/cm^ (400 psi) , -0.001 < m* < 0.031. 

Figure 8.- Heat transfer distribution over conical model 1. 
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(b) = 276 N/cm^ (400 psi), -0.001 < m* < 0.248. 

Figure 9.- Continued. 


51 







0 12 3 


B 

(a) Data generated by boundary-layer code of reference 5 for 

hemispherical model. 

Figure 10.- Correlation of laminar heating rate ratio with mass addition 

parameter 
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Figure 11.- Effect of mass addition on stagnation point heat transfer 



EXP. Re, 


ooQ 

1.2x10® 

2.4x10® 

4.8x10® 


THEORY EQ. (4) 






trans 


1 


I 

1 — 



Mw 


Figure 13.- Correlation of transition Reynolds numbers, hemispherical model 
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Figure 14.- Photomicrograph of hemispherical model surface (20x) 
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Figure 15.- PANT boundary-layer transition correlation, hemispherical model 
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CASE I. FULLY TURBULENT HEATING REFERENCE 
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CASE II. MODIFIED TURBULENT HEATING REFERENCE 


(a) Experimental data, hemispherical models. 

Figure 16.- Correlation of turbulent heating rate ratio _<^turb 

addition parameter B^urb* 


with mass 
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